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Abstract 

In this paper we investigate expanding Bianchi type I models with two tilted fluids with the same 
linear equation of state, characterized by the equation of state parameter w. Individually the fluids 
have non-zero energy fluxes w.r.t. the symmetry surfaces, but these cancel each other because of the 
Codazzi constraint. We prove that when w = the model isotropizes to the future. Using numerical 
simulations and a linear analysis we also find the asymptotic states of models with w > 0. We find 
that future isotropization occurs if and only if w < |. The results are compared to similar models 
investigated previously where the two fluids have different equation of state parameters. 
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1 Introduction 

There have been numerous investigations of spatially homogeneous (SH) cosmological models with a 
perfect fluid as a matter source. In the study of these models one usually distinguish between two 
different scenarios: when the fluid congruence is normal to the hypersurfaces of homogeneity and when 
it is not. The models of the first kind are called orthogonal or non-tilted, and the second kind are called 
tilted models. The orthogonal models have been studied extensively over the past few decades (see [I] 
and references therein), while the tilted have received detailed investigations only more recently [2j-[TTj. 
In investigations of tilted Bianchi models, the models of type I have been ignored since they are not 
compatible with any net energy flux with respect to the homogeneous hypersurfaces, something a tilted 
fluid will induce. However, it is possible to consider models of type I with tilted fluids if one have two 
or more fluids; the fluids can then all be tilted and induce energy flux individually, but the total energy 
flux must vanish. This was done for two tilted perfect fluids with linear equations of state in [12J . The 
fluids had equations of state pm = wm pm and p( 2 ) = W( 2 ) P(2), where and p^) are the energy 
densities and pressures of the respective fluids in their own rest frames, and where it was assumed that 
< ui(2) < W(i) < 1. The case = w^) was left due to lack of space and because of its qualitatively 
different behavior - this will be the topic of the present paper. 

As an example of a model with two different perfect fluids with the same equation of state one can 
imagine a universe filled with two kinds of dust (the pressure of both fluids is zero) . This is a reasonable 
model for both the baryonic and cold dark matter components of the total matter content of the universe 
today according to observations. The study of Bianchi type I models with two tilted fluids with the 
same equation of state can be viewed as the simplest anisotropic tilted model generalizing the standard 
isotropic cosmological scenario close to flatness. 

The SH Bianchi models admit a three dimensional group of isometries acting simply transitively on 
the spacelike hypersurfaces that define the surfaces of homogeneity. The models are classified by the 
Lie algebra of the isometry group, which results in a hierarchy of models of different complexity where 
the Bianchi type I model can be found on the bottom, its algebra being Abelian and obtainable from 
all the other Bianchi models by Lie algebra contractions pQ. When the Einstein field equations are 
formulated as a dynamical system the Bianchi type I models appear as a boundary of a larger state space 
describing the other Bianchi models. Describing the dynamics of this boundary is an important first step 
to understanding the dynamics of the more general models. 

The Bianchi models are sometimes interesting even when studying general inhomogeneous cosmolog- 
ical models. This is the case for example in the very early universe near an initial singularity where 
horizons form and asymptotically shrink towards the singularity which result in the equations asymptot- 
ically approaching those of homogenous models, something referred to as asymptotic silence and locality 

In this paper we analyze the system of autonomous DEs derived in [12] describing the time evolution 
of the spatially homogeneous Bianchi type I models with two perfect fluids, but in contrast to [T^] we 
study the case when the two equations of state are the same. The paper is organized as follows: In section 

2 the system of equations and constraints is derived. In section 3 the future evolution of dust models is 
described, and in section 4 the evolution of models with other equation of state is investigated. Section 
5 consists of a summary of the results and a discussion of their implication. Appendix A contains a 
description of the equilibrium points of the system and their stability properties, and appendix B studies 
the dynamics on the vacuum subset. 

2 The dynamical system 

We begin by deriving the evolution equations for the tilted two- fluid Bianchi type I models. We first give 
the evolution equations of a general Bianchi model, in terms of expansion-normalized variables defined 
relative to the timelike congruence normal to the group orbits. This derivation is described in pQ, or [14], 
but for completeness we give a short description of it here. 

First one introduces a group-invariant orthonormal frame {eo, e Q }, where eo = n is the normal to 
the group orbits and a = 1,2,3. The commutation functions, ~f a b c {t) are the basic variables: 

[e b , e c ] = 7 a f, c e a . 
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The commutator functions are normally decomposed according to 



[e ,e a ] = -[H5j + <jj + eJ 1 W]e Pl (la) 
[e a ,ep] = c 1 a pe 1 = 2a [a S (3] ' 1 + e a p n 5 ~< . (lb) 

where H is the Hubble scalar, which is related to the expansion 9 of the normal congruence n according 
to H = ^9; a a p is the shear associated with n; f2 a is the Fermi rotation which describes how the spatial 
triad rotates with respect to a gyroscopically fixed so-called Fermi frame; n af3 and a a describe the Lie 
algebra of the 3-dimensional simply transitive Lie group and determine the spatial three-curvature, see 
e.g. p. 

The energy-momentum tensor is similarly decomposed, 

Tab = pn a n b + 2q {a n b ) +p(g a b + n a n b ) + ir ab , (2) 
and is described by the source terms relative to the orthonormal frame, 

{p, p, q a , ir a p} . 

To obtain a regular, dimensionless system of equations the commutation functions and the source 
terms are normalized with the hubble scalar according to 



v _ Og.p AT _ n af3 _ a a _ fl a 

~ 3tf 2 ' 3ff 2 ' ^" ~ 3tf 2 ' aP ~ 3ff 2 ■ 1 J 

We also choose a new dimensionless time coordinate r according to 

<«> 

The evolution of H is determined by the deceleration parameter q, 

H' = -(l + q)H, (5) 

where ' denotes differentiation with respect to r. Raychaudhuri's equation gives an expression for q in 
terms of the variables (J3J): 

q = 2£ 2 + ±(fi + 3P), (6) 

where S 2 = S Q/3 . The Einstein field equations and the Jacobi identities then yields the following 

system of equations: 

Evolution equations: 

^' a fj = -(2 - g)S Q/ 3 + 2e 7l5 ( Q i? 7 - 3 TZ( a p) + 3IT a( 3 , (7a) 

A' a = [q 5 J - E a " - £f / 7 W]A p , (7b) 

(N af3 )' = [q 5^ + 2E-/ tt + 2e 7 ( Q 5 R S ]N^ . (7c) 

Constraint equations: 

= l-£ 2 + ± 3 ft-fJ, (8a) 

= (36 a ~> A p + ej^ N S0 ) E" 7 - 3Q Q , (8b) 

= A p N^ a . (8c) 

where 3r R-( a p) and 3 1Z are the trace-free and scalar parts of the Hubble-normalizcd three-curvature, 
respectively, according to: 
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3 TZ {aP) = B (a0) +2e^ {a N 0)s A r) 3 TZ = -\B<* a - 6A a A a ; B a0 = 2N ai N~< - N\ N a0 , 

and where (..) denotes trace-free symmetrization of the indices, i.e. A^ a ^ = ^4( Q( 3) — ^S a0 A a A" . The 
contracted Bianchi identities yields evolution equations for the total source variables: 

n' = (2q- l)n-3P + 2A a Q a - £ Q/3 I± a/3 , (9a) 
Q' a = -[2(1 - q) 5 J + Ej + ej\ K*] Q p + (35^ A + e a ^ N S0 ) . (9b) 

These equations are a consequence of the evolution and constraint equations ([7]) and ((8]) and give no 
additional information, but are useful as auxiliary equations. 

The system of equations (|7|) and © with ^ is not fully determined, there are no evolution equations 
for the variables R a that represent the angular velocity of the spatial frame {e a }. We can freely specify 
R a in any way most convenient. Also the evolution of the source variables P and H a0 are not determined 
until we specify a source. We consider the source to be two tilted perfect fluids, only interacting with 
each other gravitationally. We can then split the energy-momentum tensor into two parts, each satisfying 
the conservation equation separately 

T ab = T$ , V«Tg = 0, (* = 1,2), (10) 



where 



^ = iP(i)+P(i))uh jl h+P(i)9 ab - (11) 

We now impose the same linear equation of state for the two fluids, i.e., p(^\ — wp^, where w = const. 
The four velocities U(^) can be written in the form 



*(<) 



= (n° + U f i) ); n a vf {) =0. (12) 



- v a{i) v a (j) 

Since we have two separate conservation equations we now have two sets of source variables 

p w> q« n«}, (i = i,2), 

where each and Q a satisfies equations (J9J). The source variables can all be expressed in terms of the 
normalized energy density and the three- velocity of the respective fluid as 



g5 ) -(i + w )(G«r 1 ^ i) % ) , 

P (i) = vM {i) + l(l-tov)Q<$vfc, 



n 



a/3 



)(<)„« 



(13a) 
(13b) 
(13c) 



where G 



1 ± wv 2 (i) and «2U 



v a(i) v (i)- One can use equations © and (fT3")l to obtain evolution 

(14a) 



equations for the source terms fl^ and u?k : 

fi' w = (2g - 1 - 3w) n w + [(3w - 1) « (i)a - S Q/3 wgj + 2A a ] Qft , 
< W = (^V 1 [(1 - ^))(3^ - 1 - «; A" D/JW ) + (1 - «;)(^ + E/ wfo) v 0(i) 
- [T,J + e Q ^ (i? 7 + 7V 7 5 «, w )] U/3W - A a v 2 {l) . 



v a (i) 



(14b) 



For the Bianchi type I models we have A a 



N, 



a/3 



0. The Codazzi constraint (|8b[) then becomes, 



Qa = Qa + — 0, which, taken in combination with (fTB"]) . forces the 3-velocities of the two fluids to 
be anti-parallel. Kincmatically the situation is similar to that of Bianchi type I with a general magnetic 
field studied in [15] , and it is therefore natural to exploit the same mathematical structures in the present 
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problem. We therefore choose the spatial triad so that one of the frame vectors is aligned with the fluid 
velocities, which we choose to be e%, i.e. = (0,0, v^). Demanding that these conditions on vfy hold 
for all times lead to the following conditions 

R x = -E 23 , R 2 = E 3 i • (15) 

This leaves P 3 undetermined, however, we still have the freedom of arbitrary rotations in the 1-2-plane, 
which we use to set 

#3 = 0. (16) 
Following [15j . we introduce the variables E + , E^, E^, Ec according to 

E + = ±(E n +E 22 ), E 3 i+iE2 3 = >/3EAe i *, E_ + ^=E 12 = (E B + i E c )e 2i * , (17) 

v3 

where E_ = (En — E 22 )/(2-\/3), which leads to 

E 2 = Y? + + Y? A + E| + Y? C . (18) 

The angular variable <j> decouples from the other equations, <f>' — — Ec, which reduces the dimension of 
the dynamical system by one. The resulting dynamical system is the following: 



Evolution equations: 



E' + = -(2 - g)E + + 3E 2 , - Q {1) v {1) - Q {2) v {2) , (19a) 

E' A = -(2-g + 3E+ + \/3E B )E A , (19b) 

E' B = -(2-g)E B + V3Ei-2\/3E 2 7 , (19c) 

E' c = -(2-g-2V3E s )E c , (19d) 

«' W = (G^y'il - ^))(3 WW - 1 + 2E+) U(i) , (19e) 

^(i) = ( 2 9 - 1 - 3tO(i))f2(i) + (3t0(j) - 1 + 2E + )Q( i )U (i ) . (19f) 



Constraint equations: 



= 1 — s — n (1) - n (2 ) , (20a) 

= Q ( i) + Q ( 2) , (20b) 

where 

<Z = 2E 2 + ±(f> m + 3P m ) = 2-§(O n -Pm); n m = + fi (2) , P m = P (1) + P (2) . (21) 

The assumption of non- negative energy densities, f2(j) > 0, together with (j2Tj) and (|20a|) . yields that 
3<?<2- 



2.1 The state space 

The state space consists of S = {E + , E^, Eg, Ec, Um, V( 2 ), fimj ^( 2 )} subject to the constraints (|20a[) . 
(|20b[) . Both u^) = and zj 2 ^ = 1 defines invariant subsets which means that t> 2 ^, contained in the range 
(0, 1), is bounded both from above and below. The constraint (|20a[) together with the non-negativity 
condition on the energy densities puts bounds on all the other variables and hence the state space is 
bounded. We are primarily interested in the interior state space where 

o<n (1) n (2)j 0<if 4) <l, (22) 

but asymptotically the orbits of the system may approach the boundary and hence we consider its closure, 
S, which means that we consider the set defined by E 2 < 1, < w 2 ^ < 1; < O(j) < 1, in such a way that 
the constraints (f2TJ)) are satisfied. 
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The dynamical system have the following discrete symmetries: 



C 



(«(1),V(2)) ~> -(V(1),U(2)) 



(23) 



We therefore assume without loss of generality that XU e [0,1], Xc € [0,1], «(i) € [0, 1], and «( 2 ) <G [—1,0]; 
the solutions in the other sectors of the state space are easily obtained by means of the discrete symmetries. 

In |12j one had two monotone functions which here become constants of motion. Only one of them is 
independent of the Codazzi constraint, however, and can be written as 



= k. 



(24) 



where k is a positive real constant. The submanifolds defined by (|2"4"|) foliates the state space and are 
dependent on w. Projections onto Vm x W( 2 )-space for dust and almost stiff equations of state are shown 
in Figure [TJ 




0,4- 

«(2) 
-0,6 




(a) Dust, w = 0. 



0,2 0,4 „, 0,6 

1>(1) 

(b) Almost stiff, w = 0.99. 



Figure 1: The leaves of foliation in vi^ x «(2)-space. 



The constant of motion (|24[) correlates the velocities of the fluids such that if one of them becomes 
extremely tilted or orthogonal asymptotically so must the other. Note also that it is not possible that one 
of the fluids will dominate over the other asymptotically; if the normalized energy density of one of them 
vanishes then this is also the case for the other, which follows from considering the constraints (|24| and 
(|20bjl simultaneously. This excludes a large part of the boundary of S from consideration since it may not 
be approached. When examining invariant subsets and fix points on the boundary in the following we 
will only consider those that can be approached from the interior, i.e. not those where O^-j = ^ 
v (i) = ^ v (j) , or vf z) = 1 ^ vf j} , where i^j. 



2.2 Invariant subsets 



The dynamical system (fT9| . fj20[) . admits a number of invariant subsets, conveniently divided into two 
classes: (i) 'geometric subsets', i.e., sets associated with conditions on the shear and hence the metric; 
(ii) invariant sets associated with conditions on the tilt or energy densities. We will introduce a notation 
where the kernel suggests the type of subset and where a subscript, when existent, suggests the values of 
and V(2y 



Geometric subsets 



• TW: The 'twisting' subset, characterized by = 0, XU ^ 0, which leads to that the decoupled 
(^-variable satisfies <j) = const and hence X12 oc En — X 2 2- 
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Name 


< W < | 


\<w<\ 


5 <w< 1 


5 
9 


| < to < 1 


fO 
K 00 






source/saddle 




^11 






source/saddle 






center-saddle 


n (1) n (2) 

r 00 


sink 


saddle 




does not exist 


sink 


saddle 


TW n 


saddle 


TWL t)(1)1 , (2) 


does not exist 


sink 


(TWLjjjjjp] C GSt, (1 ) U(2 j) 


saddle 


GS„ (1) „ (2) 


does not exist 


sink 


does not exist 


Gii 


saddle 


(Gn c GSd( 1 )i'(2) ) 


sink 



Table 1: A list of fix points when < w < 1. The fix points are denoted by a kernel that is related to a 
subset of which the fix point belong in combination with a subscript and sometimes also a superscript. 
The subscript indicates the fix point values of vru and Vm . The superscript of the Friedmann fix points 
Fqq indicates of the values of Sim and £1(2) while it denotes the sign of E^ in the case of the Kasner 
surface KS^ . A complete characterization of the fix points is given in Appendix \K\ here we have given 
a description in terms of their stability properties. 

• 7ZT>: The constantly rotated diagonal subset, given by XU = 0, Ec 7^ (R a — 0). This subset is 
the diagonal subset, discussed next, rotated with a constant angle around e^. 

• T>: The diagonal subset, defined by E^ = Ec = 0; E# — £_, and hence R a = 0. 

• C1ZS: The locally rotationally symmetric subset. This plane symmetric subset of the diagonal 
subset is characterized by the additional condition £# = E_ = 0. This is the simplest subset 
compatible with two tilted fluids. 

Matter subsets 

• O: The orthogonal subset for which vm = V(2) = 0. In this subset f2m oc 0,(2) and the distinction 
between the two fluids becomes artificial. The orthogonal subset describes a single orthogonal fluid. 
In general this subset is expressed in a non-Fermi frame for which E^Ec" 7^ 0, however, usually 
when dealing with this case one makes a rotation to a Fermi frame in which the shear and the 
metric are diagonal so that O belongs to T>. 

• £Tn: The double extreme tilt subset where both fluids propagate with the speed of light, = 

1 = -1!(2) => £1(1) = £1(2) = 3P(i) = 3P( 2 ). 

• /C: The vacuum subset is called the Kasner subset and is defined by £l m = 0; E 2 = 1; it describes 
the Kasner solutions, but in general in a non-Fermi propagated frame, and with vu\ as test fields. 

Other invariant subsets can be obtained by taking intersections of the ones described above. 

Apart from the above subsets there are also a number of fix points. These, and their eigenvalues, 
are given in Appendix fAl but we summarize them in Table [1] along with their stability properties. If 
£l(i) and £1(2) are zer0 then vm) and «( 2 ) act as test fields. This is the case for the Kasner fix points 
K„ which all describe the same Kasner solutions; the KS^ describe a particular Kasner solution 

(l)f (2 ) 

determined by w; the line of Friedmann fix points Fqq 1 ' (2) describe the flat Friedmann-Lemaitre solution 
for a single orthogonal perfect fluid with a linear equation of state. However, all the remaining fix points 
are physically distinct. 
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3 Future isotropization of tilted multi-fluid dust models 

Arguably the most physically relevant equation of state is the case when w = 0. This describes two perfect 
fluids with zero pressure, commonly referred to as dust, and is a good description of the actual matter 
content of the universe during the matter dominated epoch. For this case we can prove the following 
proposition 

Proposition 3.1. If w=0, then the future asymptotic state of the system 119\) is the Friedmann-Lematre 
solution. 

Proof. For w = 0we have 

pn(l-if i) )nf i) ] , = 2(2g-l) 1 (25) 

which is a monotonically increasing function since q > |. Since it is also bounded from above it must 
approach a limit value and hence we have q — ► |. This implies P m — ► 0, Q m — > 1, which in turn implies 
V(i) — * 0. The constraint (|20a[) ensures isotropization and we have reduced the system to the fix points 

FLo <1)0(2) . □ 

Observations of galaxies and galaxy clusters correlates the distribution of dark matter with the bary- 
onic matter visible in the galaxies, which implies that the velocity of dark matter is aligned with the 
velocity of visible matter [16j . The Bianchi type I model is an example of how a universe where the dust 
flows initially are non-aligned can evolve into a state where they become aligned asymptotically to the 
future. 

The linear analysis of appendix A, and the dynamics on the Kasner subset, described in Appendix 
B and numerical simulations suggests that asymptotically to the past the system approaches a doubly 
extremely tilted Kasner model described by one of the fix points in (|44a[) . as conjectured in section [4~2l 



4 Future and past dynamics of general linear equations of state 
4.1 Future dynamics 

As was shown in [12 ■ no tilted two-fluid Bianchi type I models with > 0, U/U < 1 and w > | 
isotropize to the future. The theorem does not tell us what the asymptotic state is other than that it 
is anisotropic. The conclusions we make about the future global attractors for w < | rests on the local 
stability analysis of the fix points and numerical simulations, which makes us confident of the following 
conjectures: 

Conjecture 4.1. The w-limit for all orbits that have S^, T.c 7^ initially is contained in the set 
FL oo U> <2> if w <\, LRSL„ (1) „ (2) if \ < w < \, TWL„ (1) „ (2) if \<w<\, GS„ (1) „ (2) if w = § , and 
Gn if f < w < 1. 

Conjecture 4.2. The co-limit for all orbits that have T.a ^ 0, Ec* = initially is contained in the set 
FLo (1)n<2) if w < i, LRSL„ (1)U(2) if | < w < \, TWL„ (1) „ (2) if \ < w < f , and TW n if § < w < 1. 

Conjecture 4.3. The uj-limit for all orbits that have = initially is contained in the set FLg <1)S1(2) 
if w < I an d in the set LRSL„ (1)1 , (2) if | < w < 1. 

The conjectures about the future attractors can conveniently be summarized by three diagrams show- 
ing the attractors for different values of w, for the general Bianchi type I set, the 'twisting' subset and 
the diagonal subset. See figure 2. 



4.2 Past dynamics 

Numerical calculations, the local analysis of the fix points in appendix A, and the dynamics on the Kasner 
subset described in Appendix B, supports the following conjecture: 

Conjecture 4.4. The a-limit for every orbit with Q(\) > 0? w m < 1j ^2) ^ * n *^ a ^J/ on ^ e general 
geometric set with S^Sp ^ is one of the fix points on the global past attractor for the Kasner 

subset K. given in equation (1441) . 
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TWL V(1)V{2) 





GS„ (1)t , {2) 


tL oo \ , 


, Gn 



1 15 
3 2 9 



(a) The general case 



TWL 



FL 



O(l) 0( 2 ) 



1 13 

3 2 5 



1 



(b) The TW subset 



(c) The TZT>, V, CKS subsets 



Figure 2: Future global attractor bifurcation diagrams for the various geometric subsets. 



All of the tilted models approach a vacuum dominated Kasner singularity, where the influence of 
the matter becomes negligible. The fluids may either become aligned with each other and the normal 
congruence to the homogeneous hypersurfaces, or anti- aligned and extremely tilted, depending on the 
value of w and which of the Kasner states that is approached. 



5 Summary and discussion 

We have studied Bianchi type I models with two tilted fluids with the same linear equation of state, 
parameterized by the equation of state parameter w, using dynamical systems methods. The absence 
of spatial curvature in these models forces the fluids to be anti-aligned with each other. The paper was 
written with the modest ambition of completing the analysis of multi-fluid Bianchi type I models in 
[T2] where it was assumed that one the fluids were stiffer than the other (i.e. there were two different 
equations of state parameters such that > W( 2 ))- As expected the models described here are in many 
ways similar to those in [12] . and do, for example, exhibit rather similar bifurcation structure to the 
future, however , there are differences. One important difference is the existence of a constant of motion 
that correlates the velocities and energy densities of the fluids. This prohibits the asymptotic states to 
become single fluid cosmologies with an extra test field as was the case in many situations when the two 
fluids had different equations of state. 

We proved that for models with w = the system isotropizes to the future and approaches a future 
asymptotic Friedmann universe where both fluids become orthogonal to the homogenous hypersurfaces. 
For models with hi > we found equilibrium points that were future stable in the full state space 
and hence local future attractors. Numerical simulations indicate that these local attractors are also 
global future attractors, in the full state space and the invariant subspaces respectively. For equations 
of state softer than radiation the future attractors are the Friedmann points, corresponding to the self 
similar, isotropic Friedmann universe. A bifurcation occurs at the radiation value w = 1/3 where the 
future asymptotic state changes from a Friedmann universe to an anisotropic, but locally rotationally 
symmetric, state where both fluids acquire a non-zero tilt. There is a different fix point for each given 
value of a constant of motion, described in section 12.11 This result is different from the case when the 
case where one fluid is stiffer than the other, described in [12] , where the stiffer fluid becomes extremely 
tilted. A second bifurcation occurs at w — 1/2 where the future asymptotic state acquires a second 
component of shear. Since the future asymptotic fix point is contained in the twisting subset, where 
R a 7^ 0, <]) = const. , the asymptotic state is described relative to a constantly rotating frame. The frame 
is tied to the fluid three- velocity and hence indicates that the fluid is rotating relative to a Fermi frame. A 
final bifurcation occurs at w = 5/9 where the future asymptotic fix point is transferred from the twisting 
subset to the general anisotropic state space where Yi C ^ 0. This implies that the decoupled variable <f> 
is linearly decreasing (<fi' = — and hence that the frame rotation vector R a is rotating in the plane 
spanned by ei and e2. Both fluid becomes extremely tilted (i.e. \v^ \ = 1). 

The past asymptotic state is for all non-self similar solutions a fix point located on the vacuum subset, 
generically either of zero or extreme tilt, depending on the asymptotic value of one of the shear variables. 
A subset of measure zero can asymptotically to the past have intermediate tilt if the shear variable 
tends to a specific value. 

One can ask wether the anti- alignment of the fluids is a physically realistic scenario. Is there some 
mechanism which asymptotically would produce such a state from more general models? For past asymp- 
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totic behavior one could argue that there might be in many cases, since Bianchi type I models seem to 
be a part of the past attractor for very general cosmological models and for these models the fluids tend 
to become anti-aligned, but to really investigate wether this actually happens one must study two-fluid 
models in a more general context, one that at least contain Bianchi type II models. 



A Fix points and local stability analysis 



Kasner fix points: There are two circles of Kasner points and two surfaces of fix points when < w. 
The Kasner circles are characterized by E + = £ + , E# = £_, T,a = Ec = 0, £2(i) = 0(2) = 0, where £± 
are constants that satisfy £ + + E 2 . = 1, and the following values of v^y. 

K oo : v (i) = v (2) = , K ?i : v (i) = -«(2) = 1 • (26a) 

The eigenvalues for the two cases are: 

Kg,: 0; A Sa ; A Sc ; A° (1) ; A° (a) ; 3(1 - w) ; 3(1 - w) , (27a) 

K°: 0; A s , ; A Sc ; Aj (l) ; A* (3) ; 2(1 + E+), (27b) 

where 

A Sa = -(3E+ + V3E_) , A Sc =2V3S_, (28a) 

A° (i) = 3w - 1 + 2£ + , Ai w = -2(3w - 1 + 2S + )/(l - tu) . (28b) 

In the Kg, case the Codazzi constraint (|20bl) is singular and hence it cannot be locally solved; in the 
other case (|20bl) has been used to eliminate £2(i)- The zero eigenvalue corresponds to that one has a 
one-parameter set of fixed points. The two surfaces of Kasner fix points are characterized by 



KS± , : Sm - £ c = = n (2) = 0, £+ = §(1-3™), £ B =±Jl-£ 2 + , (29) 



(1-w) 





(r 











where the superscript denotes the sign of Eb , and A; is the constant defined in . The relation between 
U(i), f(2) and k constrains the three free parameters and thus gives a surface of fix points. After eliminating 
nm locally by means of the Codazzi constraint (|20bj) . the eigenvalues for the Kasner surfaces are: 

KS± i)W(2) : ; ; ; A^ ; A Sc ; 3(1 - to) , (30a) 

where again As A = — (3E+ + v3£-) > Ae c = 2v3£-, where E + , E_ take the fix point values for the 
relevant line of fix points. Here two zero eigenvalues corresponds to that one has a surface of fix points 
while the third is associated with the existence of a one parameter set of solutions that are anti-parallel 
w.r.t. each other on each side of the surface of fix points. 



Friedmann fix points: On the Friedmann subset there exists one line of fix points parameterized by 
the constant of motion k (|24[) : 

FLo (1 ' n<2) : E+ - E A = E B = E c - 0, £2 ro - 1 v (1) = , v (2) = , = Vk , (31a) 

where the superscript refers to the values of £2m) and £2(2). The associated eigenvalues are: 

FLo (1,a<2) : A 1 , 2 ,3,4 = -|(l-^); 3io-l; 0, (32a) 

were we have used the Codazzi constraint (|20b[) to eliminate the variable f( 2 ) • Two of the eigenvalues of 
Ai,2,3,4 refer to X^ A and As c . 
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Fix points on C1ZS: When \ <w there is an additional line of fix points, LRSL.„ (1)1 , (2) , which enter the 
physical state space when w = |, and move into the £7?.iS-subset with increasing values of w. The line 
intersects the foliation determined by (j2"4"|) and can be parametrised by k. We have here chosen to use 
v (i) as a parameter instead for reasons of computational simplicity. In the stiff perfect fluid limit (w = 1) 
the line merge with the coalesced Kasner surfaces. The line of fix points is characterized by: 

LRSL^^^^ . 

v v v Ifi t\ 310-1 3W-1 

E A = E s = E c = , E + = -±{3w- 1) , » (2 )%) = -- — — , - — — <«(i)<l, 

z v ' v ' 5w + 1 5w + 1 

3 (1 - io)(5u) + l)(3io - 1)(1 + CTm) 3 (1 - iu)[(5tu + l) 2 v 2 1} + w(3w - l) 2 ] 

(l + W )[(5w + l)« 2 1) + (3u;-l)] ' ^ (2) = 4 (l + w)[(5w + l)w 2 1) + (3u;-l)] ' (33a) 

After eliminating fi^) locally the eigenvalues for the LRS-line are: 



LRSL V(1)t)(2) : \ Sa =3(2^-1); A Sb = A Sc = -|(1 - w) ; 0; - to) (l ± ^A(w, , 

(34a) 

where Re Vm) < 1; since the expression for A(tU(j)) is rather messy we will refrain from giving it. 

Fix points on TW: 

TW n : E+ = -|, E c = 0, E A = E B = ^, = 1 , « {2) = -1 , = fi (a) = & . (35) 

Local elimination of Jim by means of the Codazzi constraint (120bl) yields the eigenvalues: 

6(3 - 5w) _ 6(3_5w) 
5(1 -w) ' 5(1 -w) 



A* c = i; -f; -&(i±<V5§); ^Sr^; (36) 



When ~ < u> < | there exists one more line of fix points on TW: TWL„ „ ( • This line comes into 
existence when the line LRSL 1)(1)t , (2) bifurcate into two at w = 5; it then wanders away from T> when w 
increases and eventually leaves the physical state space through TWn when w = |. Like LRSL„ C1) ^, 2 , it 
also can be parameterized by k but we choose vry\ here also for simplicity. The fix points are characterized 

by 



TWL„ (1)l)(2) 



E + = -i(3w (1) -l), E A = v /f(l- w )(2 w -l), E s = V3(2w-l) 

(l-w)(15u;-7) (l- w )(15u,-7) 
E c = 0, W(2) = __ I ____, _____<„ (1) <i 

fyi) = B(w, «(!)) , n (2) =l-J(3«;-l)(15u;-7)-fl(t«,w ( ij), (37a) 



where 



3 (1 - to) (7 - 15w)(25w 2 - 18w + 1)(1 + wvl } ) 
B ( w > v d)) = 4 (1 + W )[(_25 W 2 + 18w - l)v 2 {1) + (1 - w)(7 - 15™)] ' (3?b) 
Local elimination of On ) yields the following eigenvalues: 

TWL ni) „ (2) : A So = -|(5- 9w) ; 0; A 3 (iu, ; A 4 (w, «(i)) ; A 5 (w, ; A 6 (w, v (1) ) , 

(38a) 

where A3,4 5,6 exhibit extremely messy expressions, which we therefore refrain from giving. They all have 
the property that the real part of the eigenvalues always are negative in the domain of definition of the 
fix point set, thus making the entire line a local attractor in the range 1/2 < w < 5/9. 

Fix point in the generic geometric manifold: There exists one fix point Gn for which all the 
off-diagonal components of the shear are non-zero. It thus exists on the generic 'geometric' manifold, but 
on the 'matter boundary' £ 7~n where both fluids are extremely tilted. It is characterized by: 

G11 : E+ = -| , E A = ^J- , E B = E c = ^= , v {1) = 1 , v [2) = -1 , = fl (2) = ± . (39) 
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Local elimination of Q(t) yields the eigenvalues: 



-Vllvs. • = -I ( i ± *V 23 ± 12 ^] ; A 5, 6 = - 2 ^_^ • ( 4 o) 



At w = I there exists a triangular surface of fix points, GS^.,,,^ , connecting TWL„ (1 . t(2 . with G\\. 

GS "(i)i'(2) is s iven b y ; 



n/2 /l7w(i)W(2) + 3 v _ l 1 / 13u (i) t; (2 



1 3^3 V 4w (i) u (2) ~ 3 ' 375 ' "° 3V3\/ 4u (1) « (2) -3 

-1<«(1)«(2) <-TO' 

%) = i -^( fl + H)) 1 « (1 )(9 + 5tfo) 

3(u ( i) -U(2))(3-4w (1) i; (2) ) 3 (u (1) - u (2) )(3 - 4u (1) u (2 )J 

Local elimination of fin) yields two zero eigenvalues and four others with complicated dependence on 
and i>( 2 ) but which all have negative real part on the entire set GS^j^,,,. 

The local stability analysis of the fix points can be summarized in a table showing how the local 
attractor is transferred from set to set with increasing values of w - from the isotropic Friedmann solutions 
for sub-radiation equation of states to the increasingly anisotropic solutions when the fluid becomes stiffer. 

invariant 

subset: TCO CRS TW GENERIC MANIFOLD 

local 

sink: FL q ' ( ' — > LRSL„ (1)1J(2) > TWL^^^^j — > GS„ (1)l)(2) — > G\\ 

g[0,|) e(|,i] G(i |) § e(§,i) 

B The /C subset 

We here discuss the Kasner subset K, with the state space K = {£+, Y, A , S^, £cj um, «( 2 )}, subjected 
to the constraints (|20a[) and (I24[) . The equations for the test fields e [0,1], v^ 2 ) £ [ — 1, 0] decouple 
from those of the shear but are still coupled to each other through (pM|) . The state space therefore can 
be written as the following Cartesian product: 

K = KP x f (2) }, KP = {S + ,E i! S B ,E c }, (42) 

where KP is the projected Kasner state space, which of course is subjected to S 2 = 1. By determining 
the a- and cj-limits for solutions on KP one can then determine the asymptotic states of and i>( 2 ), 
and thus the a- and w-limits for solutions on IC. Let us therefore first turn to the equations on KP: 

E' + =3T, 2 A ; £^ = -(3£ + + V3£b)£ a ; T! b = V3Z\ - 2^31]^ ; E' c = 2V3£ s £c • (43) 

This system is defined on the compact space £ 2 = 1. Since S + is monotonically increasing on a compact 
space it must approach a constant value, hence we have T? A cx — * asymptotically both to to future 
and to the past. Since we also have (S+ — v^Es)' = — 2v / 3£ 2 7 < 0, Sc will by the same argument also 
vanish asymptotically to both the future and to the past. £,4 = £c = defines a circle of fix points, 



the projected Kasner circle: KP , see Figure 3(a) It is described by £ + = = £_, where the 

constants £+, £_ satisfy £+ + Y?_ = 1. 

From this we conclude that the the a-limits for all solutions with T, A T,c / on K.T resides on the 
local source of KP°, yielding a segment on KP° characterized by — 1 < £+ = £+ < — < £_ < 
i.e., the segment consists of sector (213) together with the fix points Q 2 and T3 on KP°. The w-limit 
resides on the local sink, which consists of segment (312) together with the fix points Q3 and T 2 , see 
Figure [3(b)1 

The a-limits for solutions on K. are determined by the a-limits on 1CP which determine the asymptotic 
limits for v {i) . The equation for € [0, 1] on KP° is given by: = (G^V 1 ! 1 - v fi))(3w - 1 + 
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(a) Kasner sectors (b) The a- and u;-limits on K.V 



Figure 3: The projected Kasner circle KP° is divided into sectors (i,j,k), defined by Ej < Ej < Efc, 
where k is a permutation of 1, 2, 3, and where Si = E + + \/3E_, E 2 = E + — V3E_, E 3 = — 2E + , and 
the points Q Q , corresponding to the non-flat plane symmetric Kasner solution, and T Q , corresponding to 
the Taub form for the Minkowski spacetime. The global past attractor for the general geometric set with 
S^E^ ^ on JCV consists of sector (213) together with Q 2 and T 3 on KP°. The global future attractor 
consists of sector (312) together with Q 3 and T 2 on KP°. 

2E + ) |u(j)|. It follows that the a-limits for all orbits on JC on the general geometric set with E^Ec ^ 
resides on the global past attractor where the subscript denotes the range of values of w, given by 



^IjKKfpS+el-i,-!]}. ( 44a ) 

A {w= * } =i K ?i : ^+ e "I)} U { KS W> : ^+ = -3> > ( 44b ) 
-4 { | <tu} ={K?i : S+ G [-1, -|(3^ - 1))} U {KS+ 1)t)(a) : S + = -I(3t« - 1)}U 

{K°:E + e(-±(3u; -!),-§]}. (44c) 



The local stability analysis of the Kasner circle in the full Bianchi type I state space shows that the entire 
circle has unstable modes and will therefore not attract any orbits outside of the Kasner subset; the sinks 
on JC are saddle points in the full state space. 
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